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Abstract. Undoped and strongly photoexcited semiconductor superlattices with 
field-dependent recombination behave as excitable or oscillatory media with spatially 
discrete nonlinear convection and diffusion. Infinitely long, dc-current-biased 
superlattices behaving as excitable media exhibit wave fronts with increasing or 
decreasing profiles, whose velocities can be calculated by means of asymptotic methods. 
These superlattices can also support pulses of the electric field. Pulses moving 
downstream with the flux of electrons can be constructed from their component 
wave fronts, whereas pulses advancing upstream do so slowly and experience saltatory 
motion: they change slowly in long intervals of time separated by fast transitions during 
which the pulses jump to the previous superlattice period. Photoexcited superlattices 
can also behave as oscillatory media and exhibit wave trains. 
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1. Introduction 

Vertical nonlinear charge transport in undoped, weakly coupled, photoexcited 
superlattices (SLs) has many similarities to extended excitable or oscillatory media. 
In both cases, traveling waves such as wave fronts, pulses, or wave trains play an 
important role in the dynamics of the system. |1] In undoped SLs, the formation 
and dynamics of electric field domains (EFDs) are the basic mechanisms behind 
more complex phenomena such as self-sustained oscillations of the current (SSOC) 
through voltage-biased SLs, chaos, etc. Interfaces separating EFDs are wave fronts 
or wave front-bound pulses, and pulses are the building blocks of wave trains. While 
there are at least two stable EFDs with an almost constant value of the electric 
field in the case of doped SLs, there may be only one stable EFD in the case of 
an undoped SL under high photoexcitation. [2] This is similar to what happens in 
excitable media such as those described by the FitzHugh-Nagumo (FHN) model for 
nerve conduction. [3l 0] For background information on this model and related models, 
see Refs. O EJ [TJ El El [TOl [TTl [1] . The case of an oscillatory medium in which there is a 
single unstable spatially homogeneous state can also be realized in both the FHN model 
and a strongly photoexcited undoped SL.[2] However, the electric field in a voltage- 
biased SL with a finite number of periods is constrained to have a fixed area underneath 
its spatial profile, and this constraint causes it to exhibit a more complex behavior. [2] 

The key ingredient to the complex dynamics of electron transport in photoexcited 
undoped SLs is the introduction of an electron-hole recombination coefficient, which 
is a function of the electric field. Nonlinear charge transport in weakly coupled 
undoped type-I superlattices (SL) under photoexcitation is well described by spatially 
discrete drift- diffusion systems (DDDS) of equations. |12l [131 E] earlier models, the 
electron-hole recombination was considered to be a constant, independent of the electric 
field. p!2] In this case, there are three stable states with spatially homogeneous field 
profiles, and the predicted nonlinear phenomena are then quite similar to the ones 
observed in the much better known case of doped SLs (in which the hole density is 
replaced by a constant doping density): the dynamics is organized by stable wave fronts 
that join the stable homogeneous states. |l4j Introducing a field-dependent electron- 
hole recombination coefficient in the model has striking consequences. [2] For high 
photoexcitation densities, it is possible to find only one stable EFD, not two as in 
the case of a constant recombination coefficient. [12 j Then, under a dc voltage bias, 
periodic or chaotic SSOC may appear. The field profile during self-oscillations may 
exhibit nucleation of dipole waves inside the sample, splitting of one wave in two, and 
motion of the resulting waves in opposite directions. Despite the strong asymmetry 
induced by the drift terms in the equations, these dipole waves resemble the pulses in 
the FHN model for nerve conduction |9l HI] and are quite different from field profiles for 
a field-independent recombination coefficient. |12j 

In this paper, we construct relevant stable solutions of the DDDS of equations 
for dc-current-biased SL by means of asymptotic and numerical methods that extend 
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and refine those used for tlie discrete FHN system. [13 [15] For an undoped SL tliat 
acts as an excitable medium (only one stable EFD), we find that electric field pulses 
moving downstream with the electron flux can be described using matched asymptotic 
expansions based on separating the sharp leading and trailing wave fronts of the pulse 
from smoother regions outside them. Pulses moving upstream do so much more slowly 
and experience a saltatory motion, in which intervals of a slow change are separated by 
fast changes during which the pulse jumps to the previous SL period. For an undoped 
SL acting as an oscillatory medium (a single unstable EFD), we construct periodic wave 
trains consisting of an infinite succession of pulses. 



2. Model equations 

Written in nondimensional form, the equations governing nonlinear charge transport in 
a weakly coupled, undoped, photoexcited GaAs/Ala;Gai_a;As SL are:[2j 

Fi - Fi_i = u{ni- Pi) , (1) 
( 

dp, 



dF- 

S + niv{F,) - D{F,){n,+i - n,) = J , (2) 



l-r{Fi)niPi, (3) 

where i takes on integer values corresponding to the spatial periods of the SL. —Fi, 
rii and Pi are the average electric field, electron and hole surface densities at the i^^ 
SL period, respectively, and the dimensionless parameters u, S, and r are described 
below. Equation ([T]) refers to the averaged Poisson equation, while Eq. ([2]) corresponds 
to Ampere's law: the total current density J equals the sum of the displacement current 
density 6dFi/dt and the electron tunneling current density across the i^^ barrier that 
separates wells i and i + I, i.e., Jj^j+i = riiV^Fi) — D{Fi){ni^i — Uj). Note that the 
tunneling current density has a discrete drift term with electron velocity v{F) and 
a diffusive term with diffusion coefficient D{F). Charge continuity is obtained by 
differentiating Eq. ([1]) with respect to time and using Eq. ([2]) in the result: 

.5|(n.-|>.) + J.^.+i-J.-i^. = 0. (4) 

Tunneling of holes is neglected so that only photogeneration and recombination with 
electrons enter Eq. ([3]). The field-dependent quantities f, D, and the recombination 
coefficient r are described in Ref. [2] and their field dependencies are depicted in 
Fig. [U The dimensionless parameter i/ is a ratio between the carrier density scale 
determined by photogeneration and recombination and the carrier density determined 
by scattering processes. The dimensionless parameter 5 is the ratio between the time 
scales of dielectric relaxation and recombination. |2] These parameters depend on the 
photoexcitation density and the Al content x G (0, 1] in the SL barriers. For high 
photoexcitation densities, 5^1, whereas 5^1 for small photoexcitation densities. 
The parameter v can be of any order. 
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l-r{F)p^ 



(5) 



It is interesting to depict the phase plane corresponding to spatially uniform 
solutions of Eqs. ([I])-(l3]) with Ui = pi = p, Fi = F: 

8'^— = J — pv{F), and 

Generically and for fixed J, the nullclines v{F)p = J and r{F)p'^ = 1 intersect in one 
or three fixed points, depending on the Al content x in the barriers, a parameter that 
controls the functions v{F), D{F) and r(F).[2] At these fixed points, 

v{F) 



j{F) = J, j{F) 



(6) 



The function j{F) is depicted in Fig. [2l For 0.45 < x < 1, there are three fixed 
points of the system in Eq. ([5]), one on each of the three branches of p = J/v{F), and the 




Figure 1. (Color online) Electron velocity v, diffusion coefficient D, and 
recombination coefficient r vs electric field F for a GaAs/Alo.sGao.yAs SL with a 
well width Lw = 10 nm and barrier width Lb ~ 4 nni. The lattice temperature is 
assumed to be 200 K. 
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Figure 2. (Color online) Local current density j vs electric field F. For a fixed value 
of the total current density J, there may be one or three zeros of j{F) — J depending 
on the Al content x. The dashed lines are the boundary condition J = aF for the 
different values of the dimensionless contact conductivity a indicated in the figure. 
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ratio between jmax — jmin and the average current (jmax + jmin)/2 is sufficiently large. In 
this case, nonlinear phenomena are quite similar to those observed in doped SL: static 
electric field domains with domain walls joining the stable branches oi p = J/v{F), 
SSOC in voltage-biased SL due to the recycling of pulses formed by two moving domain 
walls bounding a high field region, etc. For < x < 0.45, j{F) is either increasing 
for positive F (if < a; < 0.25) or the ratio 2(jmax - jmin)/(jmax + jmin) is small (if 
0.25 < X < 0.45). For < x < 0.25, there is a unique fixed point at F = F^,, which for 
an appropriate value of J may be located at any of the three branches of p = J/v{F). If 
the fixed point is located at one of the two stable branches ofp = J/v{F) for which v{F) 
has a positive slope, the dynamical system of Eq. ([5]) is excitable, whereas it is oscillatory 
if the fixed point is located on the second branch of v{F) with negative slope. We shall 
see later that quite unusual phenomena are found for these cases. Under dc current 
bias, it is possible to have pulses moving to the right or to the left and periodic wave 
trains. Under dc voltage bias, these pulses and wave trains may give rise to SSOC. [2] 



3. Wave fronts in a dc-current-biased photoexcited SL behaving as an 
excitable medium 

Wave fronts and pulses of the electric field and carrier densities are key elements in 
the description of stable solutions of our model equations. We will start by describing 
these solutions for an infinite SL under constant current bias J in the limit of high 
photoexcitation densities, 5 <^ 1, in which the dynamical behavior of the SL is richer. [2] 
In this Section, we focus our attention on wave fronts, while we will deal with pulses in 
the next Section. 



3.1. Leading order construction of wave fronts 

A wave front is a moving interface separating regions of smooth field variation on the 
time scale t. Inside the front, Fi varies rapidly on the time scale t/6. Let us eliminate 
the electron density by using Eq. ([T]) in Eqs. (l2])-(j3]): 

^-7r + n^i) D{Fi) 

at u V 

= J-v{F,)pi + D{Fi){pi+,-pi), (7) 
% =l-r{F,)p^(p^ + ^^^^^] , (8) 



at \ ^ , 

In the regions where the field varies smoothly, we can set 6 = 0, Fi = -Fj_i, Ui = pi in 
Eqs. (I7])-(l8]), thereby obtaining the reduced problem 

Piv{Fi) = J, (9) 
l-r(i^,)p?. (10) 



dt 

In the wave fronts that separate these regions, the electric field and the hole density 
vary rapidly as Fi{t) = F{C,), pi{t) = p(^), with ^ = i — ct/5. In these regions, Eqs. ([7]) 
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and ([H]) yield to leading order 



dF 



p 



V 



vim) 



^ omfJ^±3±^^Sl^2m , (n, 

Thus p is a constant equal to the value Piit) at the last point in the region of 
smooth variation before the front. Let F^^\p) < F^'^\p) < F^^\p) be the solutions 
of J/v{F) = p for Vrain < J/p < t'max, where the local maximum and minimum of 
the velocity v{F) are reached at (-Fmax, ^^max) and at (-Fmin, I'min), respectively (with 
-Fmax < -Fmin)- Let US assumc that the front velocity is positive, c > 0. Equation ( |TT|1 
has decreasing front solutions (DFs) such that the profile F{i — ct/6) is a decreasing 
function satisfying the boundary conditions: 

F i-oo) = F,{+oo) = F^'\p), F(+oo) = F,(-cx)) =F«(P) > (13) 

and increasing front solutions (IFs) with an increasing profile F{i — ct/6) that satisfies 
the boundary conditions: 

F(-oo) = F,{+oo) = F^^\p), F(+oo) = F,{-oo) = F'^^\p) . (14) 

In Fig. [HI we show the nuUclines p = J/v{F) and p = [r(F)]^^/^ for fixed 
J corresponding to Eq. ([5]) in a GaAs/AlAs SL. Depending on the parameters, the 
increasing function p = [r(F)]"^/^ may intersect the cubic p = J/v(F) in one, two or 
three critical points. At least one critical point on the first or third branch ofp = J/v{F) 
makes our system behave as an excitable medium, whereas if there is only one critical 
point on the second branch of p = J/v{F), our system has an oscillatory character, as 
we will see later. In Fig. |3l we have also shown a horizontal line for a given constant 
value of p, which joins the first and third branches of p = J/v{F). This line corresponds 
to a wave front with fixed p as explained above. 




Figure 3. (Color online) p-F phase plane exhibiting nullclines for fixed J and a wave 
front at constant p in a GaAs/AlAs SL. 



Wave fronts, pulses and wave trains in photoexcited superlattices behaving as excitable or oscillatory medial 



0.25 



0.2 



0.15 

o 

0.1 



0.05 



0.995 1 1.005 1.01 1.015 

J 

X 10"^ 



15 



10 



o 5 







-5 



1 1.005 1.01 

J 

Figure 4. (Color online) Diniensionless velocities of DF and IF as a function of 
dimensionless current density J for (a) v — 8.45 and (b) v — 330 with x = 0.25 and p 
fixed. The DF (resp. IF) velocity increases (resp. decreases) with increasing J. 





For sufficiently large i/, there are two critical values of the current density, Jd and 
Jc2, so that IFs move to the right (c > 0) if J < Jd, are pinned (c = 0) if Jd < J < Jc2, 
and move to the left (c < 0) if J > Jc2, as shown in Fig. Hl^a). A similar picture holds for 
DFs [cf. Fig. m^b)]. Near the critical currents Jd [i = 1,2), the wave front profiles and 
their velocities can be approximately found by means of the theory of active quantum 
wells (QWs) developed for doped SLs in Ref. [17] (the only change is replacing J/p and 
u/p instead of J and u in the expressions for doped SLs). In the limit u — 0+, the 
equations for the fronts can be approximated by their continuum limit, which is a ffist- 
order hyperbolic partial differential equation having shock waves among its solutions, 
and the shock velocity gives a good approximation of the wave front velocity. [IT] 

The theory we have just sketched holds in the limit as (5 — )■ 0, and we would expect 
that it also holds for sufficiently small positive values of 5. However, a comparison of the 
front velocity given by the asymptotic theory with the one obtained by direct numerical 
simulation shows a remarkable difference even for quite small values of 6. The difference 
between the approximate and numerical velocities for 6 = 0.018 is about 0.06, which 
gives a relative error of 25%. This error is reduced to 2% for 6 = 0.001. Nevertheless, we 
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would like to have an asymptotic theory, which is better than this. How do we correct 
the simple wave front construction given above? 



3.2. Corrected asymptotic theory of wave fronts 

To correct the previous leading order theory of wave fronts, we use the fact that the DF 
or IF profiles are monotone functions of a variable ^ = i — ct/5. Thus 

with < (resp. > 0) for DF (resp. IF). In either case, we can find 

^ = E{F) which solves = F . 

Then the finite differences 

- = J-(0 - -l) = F,- T{E{F,) - 1), and 
Fi+i + - 2Fi = + 1) + TiE{F,) - 1) - 2F, 

can be considered to be functions of Fj. Therefore, we can derive from the equation for 



(15) 
(16) 

(17) 



dpi 
'd^ 



1 - r{Fi)pi Pi + 



and from the equation for the wave front field profile the following equation for pi as a 
function of the field F,-: 



dpi 
dF, 



l-r(Fi)p, + 



J 



F)-F^\ v(F,) D{F,) F.+i+F»-i-2F, 
Pi H ^ ) — + — [Pt+1 - Pi H ^ 



(19) 



Now we integrate this equation and iterate the result starting from the value Pj(Fo) = p, 
thereby obtaining 



Pi^ P 



6 

+ 7 



r(F)p (p 



V 



dF 



I - (^p+ ^--^(S(F)-l) j v(F) ^ Dip ^ F{SiF)+l)+T(S{F)~l)-2F ^ 



(20) 



up to terms of order 5^. 

The starting point Fq should be selected so as to ensure convergence of the integral 
in Eq. (120|) . li p = p^: corresponds exactly to a critical point, we can select Fq as the field 
value of the same critical point F^. Then the integrand has a finite limit as F — )■ Fq, 
and the integral in Eq. (|20l) converges. If this is not the case, we need a nonzero value 
of the finite difference Fq — J-'(S(Fo) — 1) to obtain a nonzero value of the denominator 
in Eq. f|T9l) . We use the slow scale equation to calculate Fj = $(pj). Inserting this 
function in Eq. ( ITOl) . we obtain 

dpi _ 
dt 



r(<l>(pi))p2 . 



(21) 
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Figure 5. (Color online) Comparison of the numerically obtained and asymptotic 
approximations for the DF velocity as a function of ln(l/(5) using the parameter values 
X = 0.3, J = 1.006, and v = 4.449. 



We now solve this equation for an initial value Pi(0) = p to obtain pi = Pi{t\p). Then 
Fiit) = ^{pi(t;p)). Using that Fj is a function of i — at/ 6 in the wave front profile, we 
have -Fj_i(t = 0) = Fi(t = S/c) = ^{pi{—6/c;p)). The equation 



^(2(Fo) - 1: 



^ {Pi{ 



(22) 



determines Fq as a function of p with Fq — J^(E{Fq) — 1) 7^ 0. 

Once we have determined Pi as a function of Fi by means of Eq. fl20l) . we can solve 
the fast equation 

F^ 



C—+ [pm + - 



v{F.) 



D{F,)[pi+,{F,^,)-Pi{F,) + 



"i+i + Fi^i — 2Fi 



J 



(23) 



with Fi 



F^'\p), (24) 
F^'\p), (25) 



F(^), Fi±i = F{^ ± 1), and boundary conditions (for c > 0) 

Fi(-cx)) = F{+oo) = F^'^\p), and i^,(+oo) = F(-oo) 

for the IFs and 

Fi(-oo) = F(+oo) = F^^\p'), and i^,(+oo) = F(-oo) 

for the DFs. In Eq. ( 12^ . the value p' ^ p is determined by solving Eqs. ( 120|) and ( 123|) . 
until the first branch of J/v{F) is reached. Similarly, for a DF, in Eq. (!25l) is found 
by solving Eqs. ( I20l) and (l23l) . until the third branch of J/v{F) is reached. The front 
velocity is now a function of 5, and Fig. [5] shows that it is a much better approximation 
to the numerically calculated front velocity than that given by the leading order theory. 
In practice, it is easier to calculate directly the finite differences =F(-Fj — Fj-ti) as functions 
of F . 
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Figure 6. (Color online) (a) Numerically obtained field profile of a pulse moving with 
positive velocity for J = 1.013, u = 8.44565, S = 3.7288 x 10"^ (b) piy^se plane 
showing the nullclines and the motion of the 150th QW as the pulse traverses it. The 
laser intensity is 120.5 kW/cm^. 



4. Pulses moving downstream in an excitable SL 

As we see in Fig. [6]^a), a pulse moving downstream with positive velocity consists of 
regions of smooth field variation on the time scale t, separated by sharp interfaces in 
which Fi varies rapidly on the time scale t/6. To find an asymptotic approximation to 
the pulse profile, we first use the leading order description of its component wave fronts, 
according to which pi is a constant independent of i inside the wave front. We can now 
discuss different regions in the asymptotic description of a pulse, recalling that the field 
profile is the mirror image of the motion of a QW, Fi{t) = F{i — ct/6). 

4-1. Only one critical point on a stable branch of J/v{F) (leading order theory) 

First of all, we shall describe the pulses for values of J and u such that there is one 
critical point {F^,p^) located on the first branch of J/v{F). There are no other critical 
points or, if they exist, they are located on the second branch of J/v{F). Similar results 
are obtained if the critical point is on the third branch of J/v{F). 

• The region of smooth variation of F in front of the pulse is described by Eqs. ([9]) 
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and (fTUj) . In this region, Fj = F^^\pi) so that 

|^ = l-r(F«(p.))p^ 

and initial data evolve exponentially fast toward equilibrium, Fi = F^,, pi = p^,. 

• For the pulse leading edge, let p{t) be the value of pi at the last point of the region 
in front of the pulse. Eventually, p p^,. At the leading edge, Fi{t) = F{i — ct/5) is 
a DF moving toward the right with velocity C = c{J, v,p)/5 measured in QWs per 
unit time t. The DF satisfies Eq. (ITT]) with the boundary conditions in Eq. (fT3l) . It 
is convenient to call c_(p) = c(J, i/, p). Eventually, C ~ c_{p^)/6, and Fj decreases 
from Fi = Fjnax to Fi = across the leading edge of the pulse. 

• For the region between fronts, F^ = F^^\pi) and 

^ = l-r(F(^)fe)) pI 

There is a finite number of points in this region. On its far right, Pi = p ^ p^. As we 
move toward the left, pt increases until it reaches a certain value P{t) corresponding 
to that in the trailing wave front. 

• For the trailing wave front, Pi{t) = p{^) = P, and Fi{t) = F{C,) is an IF satisfying 
Eq. (fTT!) with the boundary conditions of Eq. (fTi|) . This front moves with velocity 
C = c{J,h',P)/6 measured in QWs per unit time t. It is convenient to denote 
c+(P) = c{J,v,P). We shall indicate how to determine P below. Clearly, if the 
pulse moves rigidly, we should have c+(P) = c_(p*) after a sufficiently long transient 
period. 

• For the pulse tail, we again have Fj = F^^\pi) and dpi/dt = 1 — r[F^^\pi))pl. 
Sufficiently far to the left, pi = p^,, Fi = F^,. 

The number of QWs between wave fronts of the pulse can be calculated as 
follows. [T6] Let r be the delay between fronts, i.e. the time elapsed from the moment, 
at which the leading front traverses the QW i = I, to the moment, when the trailing 
front is at ? = J. Clearly, 

"^('^ dv 

(26) 



The number of QWs between fronts, m(t), is 

m = - / c^{p{t))dt. (27) 

Jt~r 

However, the separation between fronts satisfies the equation 

dm c.{p{t)) - c+{P{t)) 

H = 6 ' • ^^^^ 

The three Eqs. fl2^ . flTTI) . and (1251) can be solved to obtain the three unknowns r, m 
and P{t). The function p{t) is determined by solving Eq. (ITO!) with Fj = F^^\pi) in the 
region to the left of the leading front. 
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After a transient period, p{t) — )■ p^ and P(t) — t- P (a constant value) so that we 
arrive at the simpler expressions 

dv 

^ and (29) 

(30) 



1 - r (F(3)(p)) p2 
dm C-{p^) — c+{P) 



dt S 

instead of Eqs. fj26|) and fl28l) . respectively. The number of points at the pulse top is 
now 

_ c-(p,)r ^ c,(p,) c^p 

5 54 l-r(F(3)(p))p2 • ^^^^ 

This equation yields P as a function of m. Then Eq. ( l30l) becomes an autonomous 
differential equation for m that has a stable constant solution at m = m* such 
that c_(p*) = c+(P(m)): At m = m*, the right hand side of Eq. (!30!) has a slope 
_[1 _ r(p(3)(j/p)) p2] c;(P)/c_(p,) < 0. 

Recapitulating, for appropriate initial conditions, leading and trailing fronts of a 
pulse evolve until m reaches its stable value at which c_(p=k) = c+(P(m*)) and Eq. ( !3TI) 
holds. To compute m*, we first determine P* = P(m*) by using c_(p,,) = c+(P(m*)). 
Then we calculate r = r* (which does not depend on 6\) from Eq. (!29|) and 
m* = C-{p^)t* /6. Our construction breaks down if the number of QWs between fronts 
falls below 1. This yields an upper bound for the critical value of 6 above which pulse 
propagation fails: 6c ~ c_(]9*)r*. Pulse propagation may also fail if the leading front 
becomes pinned for a current density on the interval {Jci, Jci) as mentioned in the 
previous Section. P [HI [TBI HH] 

To calculate the asymptotic length of the pulse tail, we cannot use Eq. ([29]) with F^^^ 
replacing P^^^ in that formula, because the resulting time becomes infinite. However, 
we can calculate the time it takes for a solution p{t) of Eq. ( 12T1) with $(p) = F'^^\p) 
and p(0) = P to reach a neighborhood of p^. From that equation, we obtain 



t 



p 



r{F(^){p))p'^ - 1 



r{F^^'){p))p'^ — 1 {rp'^)i{p — p*) 



(rp"^)'* \P-P 



p 



1 



with {rp^y^ = d[r{F^^\p))p'^]/dp\p=p^ > 0, and therefore 

Jp, [r(P(i)(p))p2 - 1 (^rp^y^{p - p^) _ 

for long times such that p{t) is sufficiently close to p*. The time needed for p{t) to go 
from P to + (P — p*)/e is then = T + l/(rp^)'^. The tail length is approximately 
given by M = c_{p^)T^/5. 

4-2. Only one critical point on a stable branch of J/v{F) (corrected theory) 
How do we correct this pulse construction using our improved theory of wave fronts? 
• The region of smooth variation in front of the pulse is as described above. 
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• However, we have to use Eq. (1251) instead of Eq. (ITT!) including the boundary 
conditions of Eq. fl2^ or fl25l) to construct the leading and trailing wave fronts 
of the pulse. Let us assume that the pulse moves from left to right. The leading 
wave front is a DP moving with velocity c_(p)/5 with c_(p) = c(J, z/, p, 6) given by 
Eq. (I23ll and the boundary conditions of Eq. (I25ll . While p{t) is the value at the 
initial field on the first branch of J/v{F), p'{t) is the hole density at the final QW 
of the DF which is on the third branch of J/v{F). The time it takes for a QW to 
move from (F^^\p),p) to {F^^\p'),p') is of order S, and it needs to be considered 
when constructing the pulse. 

• In the region between leading and trailing fronts, Fi = F^^^pi). On its far right. 
Pi = p' —> p'^, where we call pl^ the value of pi corresponding to p = p^,. As we move 
toward the left, Pi increases until it reaches a certain value Pit) corresponding to 
that in the trailing wave front. 

• The trailing wave front is an IF moving with velocity c+(P)/5, with C-(.(P) = 
c(J, z/, P, (5) given by Eq. ( l23l) and the boundary conditions of Eq. ( 12^ . The hole 
densities at the initial and final QWs of the IF are P and P', respectively. The 
corresponding fields are F^^\P) and F^^\P'), respectively. Again the time it takes 
for a QW to move from {F^^^{P),P) to (p(i)(P'), P') is of order 6 and will be 
included in our calculations to order 6. 

• The pulse tail is as described above except that its first QW has a hole density P' 
instead of P. 

Equations ( l26l) to ( l28l) become 

dp 

p'it-r) l-r(P(3)(p))p2 ' (33) 

1 /■* 

m = -y c_(p(t))dt, (34) 
dm _ c4p'{t))-c+iPit)) 

li ~ ~s ■ ^^^^ 

After the transient period, these equations become Eqs. ( 129|) - (!3T1) with p' (now time- 
independent) instead of p^. The rest of the considerations made above apply except 
that now r* depends on 6 because p'^ does. 

4-3. Two critical points on the stable branches of J/v{F) 

In the case illustrated in Fig. [TJ there is one critical point on each of the three branches 
of J/v{F). Let (P^,,p*) and {F*,p*) be the critical points on the first and third branches 
of J/v{F), respectively. This case is very similar to that described before for the case 
of a sole critical point on the first branch. However, now c_(p^,) ^ c^{p*) except for 
particular values of J. Therefore, these pulses do not move rigidly in general: they will 
shrink and disappear or grow indefinitely. Regions 1 to 3 of the pulse in Fig. Ul^a) are 
identical to those of the pulse with only one critical point, except that now p ^ p^ and 
P ^ p*, and m either decays to zero or it grows indefinitely. 
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Figure 7. (Color online) (a) Numerically obtained field profile of a pulse moving with 
positive velocity when there are three critical points in the phase plane, (b) Phase 
plane showing the nuUclines and the motion of the 150"* QW as the pulse traverses it. 
J = 1.009. The photoexcitation intensity is 120.5 kW/cm^. 



5. Wave trains moving downstream in a dc-current-biased photoexcited SL 

For the spatially discrete FHN system, wave trains were constructed by A. Carpio using 
matched asymptotic expansions. |T5] A similar construction could be carried out using 
our improved theory of wave fronts. A wave train consists of a periodic profile F[^), 
p{^), ^ = i — ct/S, with period L and velocity c. Figure El^a) gives the corresponding field 
profile: a first smaller pulse triggers a periodic succession of equal pulses, which become 
the wave train. Figure MJo) shows the passage of a wave train through a QW: starting 
at an unstable stationary state, trajectories in the phase plane evolve toward a stable 
limit cycle. The wave train profiles can be reconstructed from a time periodic solution 
Fi{t) = F{i — ct/5), Piif) = p{i — ct/6), with time period T = L6/c and velocity c. The 
spatial structure of the wave train at each fixed time t is given by Fi{t) = F{i — ct/6) 
and pi(t) = p{i — ct/6). The points contained in a period satisfy < i — ct/6 < L, that 
is, ct/6 < i < L + ct/6. As time grows, the discrete points travel along the continuous 
wave profile and are transferred from one period to the next. The number of integers i 
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Figure 8. (Color online) (a) Numerically obtained field profile of a pulse moving with 
positive velocity followed by a wave train when there is only one critical point in the 
phase plane, (b) Phase plane showing the nuUclines and the motion of the 350*'' QW 
as the pulse traverses it for J = 1.0165. The photoexcitation intensity corresponds to 
120.5 kW/cm^. 



one can fit in an interval of lengtli L is tlie integer part of L. 

Assuming a positive velocity, there are four stages in one time period of the wave 
trains (Fi(T, 5),pj(T, 5)): 

• For leading-edge DFs given by Eqs. ( l20l) and ( l23l) using the boundary conditions 
of Eq. (1251) . the DP velocity is C-{p)/5 with c_(p) = c{J,h',p,6), where p(t) is the 
value at the initial field on the first branch of J/v{F) and p'{t) is the hole density 
at the final QW of the DF, whose corresponding field is on the third branch of 
J/v(F). The time it takes for a QW to move from (F^^\p),p) to {F^^\p'),p') is of 
order 6, and we should consider it when constructing the train. 

• In the peak region between leading and trailing fronts, Fi = F^^\pi). On its far 
right. Pi = p' . As we move toward the left, pi increases until it reaches a certain 
value P(t) corresponding to that in the trailing wave front. Pit) will be calculated 
later. The duration of this stage is 
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Note that Tp may become infinite if tliere is a fixed point on the third 

branch of J/v{F). In such a case, we may have wave trains only if p < p*. 

• The traihng wave front is an IF moving with velocity c^{P)/6, with c+(P) = 
c{J,v, P,5) given by Eq. (!22l) and the boundary conditions of Eq. (!24l) . The hole 
densities at the initial and final QWs of the IF are P and P', respectively, and 
the corresponding fields are F^^\P) and F^^\P'), respectively. Again the time it 
takes for a QW to move from {F^^'>{P),P) to {F^^^P'), P') is of order 6 and will 
be ignored. For the wave train to move rigidly, we must have c_(p) = c+(P) = c. 
This condition gives P as a function of p. 

• In the tail of the wave train, Fi = F^^\pi). On the far right of the tail region. 
Pi = P' and Pi decreases until it reaches the value p. The duration of this stage is 



Note that Tj may become infinite if there is a fixed point {F^,p^) on the first branch 
of J/v{F). In such a case, we may have wave trains only if p > p^. 

The previous construction of the wave train gives its velocity and period T Tp + Tt 
as functions of the parameter p. The spatial period L is the integer part of cT/6. The 
number of QWs in the peak (resp. the tail) region is the integer part of c^{p)Tp/5 (resp. 
c^{p)Tt/6). Our construction clearly fails if the number of QWs in the peak region is 
smaller than one, i.e., if c_(p)Tp < 6. Thus, 6 has to be larger than c^{p)Tp. 

6. Pulses moving upstream in a dc-current-biased photoexcited SL that 
behaves as an excitable medium 

Numerical simulations of the complete model show pulses moving upstream with a 
negative velocity, both under dc current bias and under dc voltage bias. Although the 
field profile of these pulses is quite similar to that of downstream moving pulses (cf. 
Figs. [6] and [9]), there are fundamental differences between them. The velocity of these 
pulses is much smaller than that of downstream moving pulses, and they cannot be 
approximated by one IF and one DF plus regions of slow variation of the electric field. 
In fact. Fig. |4](b) shows that, contrary to the case of wave fronts moving with positive 
velocity depicted in Fig. Hl^a), it is not possible for a DF and an IF to move with the 
same negative velocity at a fixed J. Thus, our construction of pulses in Section H] cannot 
describe pulses that move rigidly upstream with negative velocity. 

The asymptotic construction of pulses moving rigidly with negative velocity under 
dc current bias is necessarily different from the case of downstream moving pulses. 
The rigid upstream motion of pulses is saltatory: although the velocity of the pulse is 
constant (c < 0), there are periods in which the QWs move slowly on the pulse field 
profile separated by fast transitions as shown in Fig. [9](a). Figure [9](b) shows that the 
hole density changes smoothly even at the fast transitions in the field profile. During 
the slow periods, the field and hole density at the QWs evolve according to the following 




(37) 
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Figure 9. (Color online) (a) Numerically obtained field profile of a pulse moving with 
negative velocity for J = 1.009, v = 70.5, 5 = 4.05 x 10~^. (b) Numerically obtained 
hole density profile of the same pulse. The laser intensity is 479.7 kW/cm^. 



equations: 

Pi H I v{Fi) - \pi+i -Pi^ 1 D{Fi) 

= J , (38) 
^ = l-r(F,)p,(p, + ^^^) , (39) 

which have been obtained by setting 5 = in Eqs. (jTl)-®- We have to solve Eq. ( l38l) for 
the profile {Fi} in terms of the instantaneous values of the {pi} and insert the result in 
Eq. (139|) . It turns out that there are several possible solutions corresponding to integer 
shifts of the pulse profile i — )• i + m, m = 0, ±1, . . .. The implicit function theorem 
establishes that, starting from an appropriate pulse-like initial condition, it is possible 
to find Fi = Fi{{pj}), provided the Jacobian determinant corresponding to Eq. fl55i) is 
not zero. This condition holds until the QWs reach the points of abrupt field change in 
Figs. M^a) or [TO](a) (phase plane) at times t = ta. Then the Jacobian vanishes, and the 
field values change according to 

^ = J+ -P. + DiF,) 
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Figure 10. (Color online) (a) Phase plane showing the nullclines and the motion of the 
8*'' QW as the pulse traverses it for the profiles depicted in Fig. HI (b) Comparison of 
the pulse velocity obtained from leading order perturbation theory and that obtained 
from numerical simulations. 



- (p^ + v{F,) , (40) 

with T = (t — ta)/S and {pi} frozen at their values at ta given by Eqs. (j38l)-(l39l). As 
explained before, as r — )• — cxd the {Fj(r)} tends to the solution of Eq. (1381) for t = ta, 
and it tends to the same profile shifted one step to the left: z— — lasr— )-+cxd. 
Then another slow stage follows. Figure [TOTb) shows the velocity of the pulse. It is 
approximately the inverse time one QW spends in the short interval of p between the 
abrupt long jump in F from the first branch of J/v{F) to a value below that in the 
third branch of J/v{F) and the abrupt short jump in F to the third branch of J/v{F), 
which occurs for a somewhat higher value of p as shown in Fig. [TOT a). 

Figures [TTT a) and [TTT b) compare the field and hole density profiles of the 
reconstructed pulse (moving with positive velocity) and those obtained by direct 
numerical solution of the equations. Similarly, Figures HS^a) and [T2r b) compare the 
field and hole density profiles of the reconstructed pulse (moving with negative velocity) 
and those obtained by direct numerical solution of the equations. The agreement is 
quite reasonable and it could be improved by finding a corrected theory of the pulse 
velocity. 
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Figure 11. (Color online) (a) Reconstruction of a pulse of the electric field moving 
with a positive velocity for J = 1.009, v = 70.5, 6 = 4.05 x 10"'' as in Fig. [71 (b) 
Reconstruction of the corresponding hole density profile. 



7. Conclusions 

An undoped, strongly photoexcited, weakly coupled, dc-current-biased SL exhibits a 
variety of pulse and wave front solutions. They can be constructed using matched 
asymptotic expansions that exploit the large separation of time scales (measured by 
the dimensionless parameter 5) between the carrier dynamics and the evolution of the 
electric field. On the time scale of the electron and hole densities, the field follows 
adiabatically the profiles thereof. 

For large photoexcitation, the field profile of a pulse typically consists of slowly 
varying regions separated by sharp wave fronts. As in the case of the FitzHugh-Nagumo 
model of nerve conduction, an asymptotic reconstruction of the downstream-moving 
pulse to leading order in 5 determines its velocity as the value for which the velocities 
of the leading and trailing fronts are the same. The comparison of this approximate 
velocity to that obtained by direct solution of the full model equations is only fair. For 
a better agreement, we have improved the description of the wave fronts comprising 
the pulse by including order 5 corrections. The resulting corrected theory describes 
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Figure 12. (Color online) (a) Reconstruction of a pulse of the electric field moving 
with a negative velocity for J = 1.009, v — 70.5, 6 — 4.05 x 10^'*. (b) Reconstruction 
of the corresponding hole density profile. 

the pulses and their velocities much better. Pulses moving upstream the electron flow 
exhibit saltatory motion: they alternate periods of very small velocity with fast motion 
during short time intervals. We have reconstructed them using matched asymptotic 
expansions in the limit of small 5, but the resulting theory is quite different from that 
of downstream-moving pulses. 
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